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TABLE OF SYMBOLS AND ABBREVIATIONS 


Definition 


Geometric buckling 


Density of the ith group of delayed 
neutron precursors 


Diffusion Coefficient for core material 


multiplication factor for infinite 
medium 


Effective multiplication factor for 


finite medium 

Reactivity 

Reactivity 

Reactivity from à disturbance 


Reactivity produced from rod movement 
due to controller action 


Thermal diffusion length of core 
materíal 


mean lifetime of thermal neutrons in an 
infinite reactor 


mean lifetime of thermal neutrons in a 
finite reactor 


Neutron density 


Normalized neutron density, 


n(t)/n (o) 


The variation of the density from its 
initial value, n(t) - n(o) 

Same as n(o), the initial steady value 
of neutron density. Also the set 
reference, the desired power level. 


Resonance escape probability 


Reactor period 





reactivity units 


reactivity units 


none 


sec 








Symbol 


Definition 


Real time 

Mean velocity of thermal neutrons 
Delayed neutron fraction 
Laplacian operator 


Average delayed neutron effectiveness 
for the six groups 


Effectiveness for the ith group 

A dimensionless parameter which appears 
in the derivation of expression for 
ker (See Appendix A) 


Fast fission factor 


Average number of neutrons emitted 
per fission (2,46 for 0235; 


Decay constant for the ith group of 
delayed neutron emitters 


Neutron flux 


Energy-averaged macroscopic neutron 
absorption cross section 


Fission cross section 
"Sigma", the summation Symbol 
Fermi age for delayed neutrons 


Fermi age for prompt neutrons 


Prime denotes analog computer voltage 
representation of quantity indicated 


none 





none 


none 


d e o 
neutrons/fission 


-1 


ses 


2 
n/em -sec, 


-1 
cm 








1. Introduction. 


The AGN-201 reactor is a low power reactor which was produced in 





quantity by Aerojet-General Nucleonics for use in education, research, 


medical diagnosis, and industrial process control [7] > The 25 by 25 cm. 


cylindrical core is constructed of discs of a homogeneous mixture of 


polyethylene moderator and UO, which is 20% enriched in vii", 

























The 


critical mass is 656 grams of y 235 dispursed in approximately 11 Kg. of 


polyethylene. Control of the reactor is accomplished by four fuel-loaded 
control rods, two act as safety rods, and two operate as control rods. 
The original licensed power was 0.1 watt but shielding modifications now 
permit the reactor at the U. S. Naval Postgraduate School to operate at 


1000 watts. At this power, a maximum flux of about 5 x 1010 


aiu cse: 
is available in a "glory hole" which penetrates the core to its center, 

The reactor is used by the Postgraduate School for a variety of 
purposes, One of the largest Gies is for the neutron activation of vari- 
ous chemical elements. This operation involves the maintenance of a 


steady power level for extended periods of time. An automatic system 


` capable of maintaining the neutron flux (reactor power) at a preset level 
would reduce operator fatigue and free him from the simple task of man- 
"ipulating a switch to keep a meter reading constant. This automatic opera- 
tion would permit him to observe more closely the ovata behavior of the 

| reactor. The objective of this thesis work was to design such a control 
system, The approach to the design was to investigate the time behavior 
“of the AGN-201, simulate that time behavior on an analog PM and use 


‘the analog simulation as an experimental tool in the design of the control 


system, 











The specifications for the automatic flux control system, here- 


after abbreviated as the AFC system, were established as follows: 


Ls 


3. 


The system will be capable of maintaining stable opera- 
tíon at any desired flux (power) level between 0.01 watts 
and 200 watts. 

Accuracy of control will be within £1% of a set power 
level. 

The neutron detection systems and control and safety 
circuits of the reactor will not be disturbed or changed 
by the AFC system, 

Control of the reactor will be effected by using the ex- 
isting control rod actuators and the rates of travel of 
the rods will not be changed. This specification stems 
from an agreement with the Atomic Energy Commission re- 
garding the proposed AFC system. | 
Override of the AFC system will be possible by manipulat- 
ing the manual dontrola of the reactor. 

The AFC system will be made as safe às possible against 
failure of its components. That ís, failure of a vacuum. 


tube or relay will not cause insertion of pesitive re- 


activity into the reactor. 











2. Type of System and Design Methods Used, 
The control problem may be explained by reference to Fig. 1. See 


Section 5 for a detailed description of the AFC system. 







REACTOR 





KINETICS 


Fig. 1 Block Diagram of the Control Problem. 
It is the function of the AFC system to compare the reactor neutron úseky“ 
n, with a reference level, no If an error is found to exist, the system 
will initiate control rod action so that rod reactivity, ke is changed 
in the proper direction to bring the error to Sepp: The possibility ex- 
ists for a reactivity disturbance, Sk, to be introduced during the 
‘performance of an experiment on the reactor, Regardless of the inter- 
action of Sk, and Sx, » the actuator must be designed to Yr the 


error as long as an error exists, 


The system specifications described above indicate the desirability 





of using a discontinuous, on-off, type control system, Specification é 


1 
Neutron density, flux and power will be used interchangeably ix 
this report. 


































arly restricting. In order to use a continuous controller 
nave the capability of varying the corrective action, As applied 
AGN-201, this means variable rod speed. Without this capability, 
ce is limited to a form of discontinuous controller which can 
a satisfactory manner using fixed control rod speed. Actual- 
| control rod on the AGN has two speeds but this is not enough 
ty for a continuous control system, A continuous controller 
ontinuous corrective action [4] + This would be unwise for a. 
od: drive mechanism, like that on the AGN-201, that was not 
lly designed for continuous duty. The AGN rod drive motors are 
- 35UbC, Shunt motors, rated for intermittent service only. In 
hg to the specification to avoid changes in the rod actuator, a 
onsideration is the fact that a power A would be requir- 
continuous control system. The addition of an amplidyne or 
ification device would require na dor and expensive changes in 
wiring of the reactor, A discontinuous controller does not re- 
ower amplifier. Amplification is effectively provided through 
on of relays. 
de from these specific considerations, the desirability of using 
nuous contactor type system for the control of nuclear reactors 
n llustrated in the past [16] [17] [23] . For the reasons des- 
ove, it was decided to design a discontinuous, on-off, type con- 
tem, 
e analytical approach to the design of a control system has been 
ed in detail in textbooks on control [9] [24] . Working with 
Esrentíal eguations or transfer functíon whích describes the be- 


of the plant, the engineer can analytically design a controller to 








make a system perform as required. Design methods favoring experi- 
mental techniques are also used in control work. After a plant trans- 
fer function has been determined or a series of differential equations 
devised that apply to the plant, a simulation devíce can be fashioned 
that will behave in a manner similar to the plant which is to be con- 
trolled. The actual equipment of the proposed control system can be 
operated with the simulation device and direct evaluation of the ability 
of the system to meet the design specification can be made. System para- 
meter changes can be readily tested with the simulation thereby simplify- 
ing the optimization of the control system design. 

In the investigation at hand, there were several factors which tend- 
ed to favor: the éxpdvinéntal approach to the design problem, The time 
behavior of a PPP reactor is nonlinear. That is, one of the seven 
differential equations which describe the kinetic behavior of a reactor 
contains a nonlinear term. In Section 3, Reactor kinetic equations, 
‘equation (2) shows the product of the multiplication factor, koff and 
"neutron density, n. Neutron density is dependent upon Kee: À trans- 
‘fer. function, by definition, describes linear behavior. If a system is 
ito be treated by analytical methods using transfer functions, it must be 
made linear. Using small signal techniques, the reactor equations can 
‘be made linear and a transfer function obtained [2] «: The result is a 
seventh order polynomial which describes the behavior of the reactor 
reasonably well under conditions of small disturbances [23] . For 
(large disturbances, the linear approximation introduces gross errors in 
'the time response. No approximations of this type are necessary in a 


plant simulation. Specifically, an electronic representation of the 


























differentíal equations of a nuclear reactor can easily handle the non- 
linear term described above by the use of a function multiplier. 

The choice of a discontinuous type controller also influences the 
design approach to be used. Describing function techniques can be used 
to approximate relay behavior in order that linear design methods may be 
applied to the problen [17] [25] > However, the method based on 
describing functions provides little information about — behavior 
ho] . hada plane analysis has been applied to the design of nonlinear 
systems but its application is limited, for computational reasons, to 
second order systems ps] . Using a Simulated control system or the 

. actual control system itself ín a closed loop with the plant simulation, 
eliminates the necessity of the above approximetions and offers the de- 

signer a direct, accurate measurement of the behavior of the actual con- 
trol system. 

Fig. 1 shows two inputs to the system, no and Šk- In the analysis 
of such a system the effect of more than one input on the system behavior 
is determined by considering the effect of each input separately and then 
applying the principle of superposition. g . This approach presumes 
prior- -application of linearization approximations to a non-linear problem, 
With anexperimental design approach, using a plant simulation, no serious 
© problems are posed by a variety of inputs. 

To avoid the approximations described above and to keep the complex- 


ity of the design problem at a minimum, it was decided to apply the símu- 


“lation approach to the design. Specifically, the differential equations 
describing the kinetic behavior of the AGN-201 were programmed on the elec- 


tronic analog computer and the computer simulation was used as a design 





tool in choosing and optimizing the components for a RET system. 

| In order to have reasonable assurance that the simulation accurate- 
ly represented the actual behavior of the AGN-201, two methods of check- 
ing the simulation were used; comparison with a digital solution of the 


reactor kinetic equations, and comparison with the actual behavior of the 





reactor itself. One might presume that comparison with the reactor would 
be sufficient proof of the validity of the simulation, Düfortunstely 
exact comparison was not possible for several reasons. One reason was 
related to the problem of applying a true step input of Fatini How- 
ever, a reasonable "forced fit", described in Section 4, was made and the 
accuracy of the analog simulation wás established. 

Utilization of a device that simulated the reactor behavior offered 
several advantages to the designer over the use of the reactor itself. 
Foremost was the lack of problems of nuclear safety. The nuclear chain 
enee requires constant supervision and control. Many control devices, ` 
te used improperly, go into a saturated condition or damage a readily re 
placable component. A nuclear reactor however, if allowed to get out of 
control, presents radiological safety hazards to SN can do 
costly damage to itself. A second advantage of a simulated reactor is 
ease of operation. An analog computer can be stopped, reset and restart- 
ed without delay. A nuclear reactor, on the other hand, requires elabor- 
ate checkout procedures prior to starting and requíres a considerable 
amount of time to shut down or return to a stable power level for a new 
run. After the initial set-up is accomplished, an analog computer can be 
operated repeatedly without any delay other than for the setting of a 


potentiometer or changing some other controller parameter for a new run. 








3. Reactor Kinetic Equations. 
The nuclear reactor kinetic equations used in this report are de- 


rived in this section. The developed equations are shown here; 
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i = 1,2, 3....6 

For definitions of the terms used, the toader is referred to the 
Table of Symbols and Abbreviations, page v. 

The above equations are simply the expressions describing the time 
(behavior of the thermal neutrons in the reactor. These expressions can 
be derived by considering the sources (production) and sinks (losses) of 
thermal neutrons. Thus, in equation form, per unit volume, 

Rate of change of neutron density - Rate of production - 

Rate of losses. 

The two sources of thermal neutrons are the prompt and delayed neu- 
-trons that have been slowed down. Also there are two SE processes 


causing loss of neutrons in a reactor, absorption and leakage’. We may 


Vibaoeptfan refers to the ica? reaction whereby a neutron is 
captured by the nucleus of some material present. Some of the captures 
^ result in fission, others produce gamma rays or other non-productive 
„results. Leakage refers to the escape of neutrons from the reactor 
System without reacting with (being captured by) a nucleus. 


















now write 


dn/dt = S *S,- LR - AR 
p d 


where 
$5 = production rate from prompt neutrons 
Sa = production rate from delayed neutrons 
LR = leakage rate of neutrons 
AR = absorption rate of neutrons 


The four terms will be described. Rather than starting from first prin- 
ciples, it will be assumed that the reader is familiar with the funda- 


mentals of reactor theory as contained in the textbook by Glasstone & 


Edlund [14] . This reference should be consulted for more detailed 


information, From reactor theory, the neutron behavior is both time and 


space dependent. The space-independent form of the expressions will be 
used, Therefore we will concern ourselves only with the time varying 


aspects of the quantities. The production rate from prompt neutrons is 


given by the expression, 
P 


- BT; 

So-(i-Ok.f.c $ 
The production rate from delayed neutrons is, 

S b -BT 

d d Pie” “AC 
As 
-B°Tp yb 
| = Pe 4 V AC 

p | AS 
: where 次 is the delayed neutron effectiveness and is defined by the ex- 


opression — B*Td, 
C 


{i= ars 一 
e P See Appendix A 


We note that P =P since the delayed neutrons are emitted at energies 


above the major resonances of Uu. 


The leakage and absorption rates for the neutrons are, 





LR — D ^p 
AR - Za 中 


The combined equation becomes 


_ ged Sat k 
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(1) 
2 
=(=DY 9) — 29 
This equation can be arranged in a more useful form by applying the 


following definitions m i 


5 l e, 
$=nv L= Ey 一 /了 人 zBz 


(this follows from space - 


D 2 
Z = Zo.” i vo = -B $ independence assumption) 
| 2 Wn B* Tp 
hou gs A 
e = / 4 L*B* See Appendix A 


Substitution of these expressions into eq. (1) and simplifying, 


dn. n Da BL) 
T d a | 


dt (2) 


2 
+ pe Yu c 
Equation (2) conteins seven unknowns, the six delayed neutron pre- 
cursors (the C's) and the neutron density (n). Six more equations are 
„obtained from the rates of change of the precursors. In equation form, 


dC, /dt - Production Rate - Decay Rate 
i = 1l, 2, ecc 
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Making substitutions similar to those made in eq. (1), 


(4) 





e tod en Ka 
at 7 r E Pe” Tp 
ic 1,2,....6 


T has a value of 1.0009 for y ?35 with y- 1.14 and will be approxi- 


mated by unity in all further expressions. 
The seven equations made up of (2) and (4) are the kinetic equations 


of the nuclear reactor. These equations will be used in Section 4 in 


setting up the analog computer solution and in programming the digital 


solution, . 
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4, Simulation of the Nuclear Reactor, 

The actual design and testing of components of the flux controller 
was accomplished using an analog -— simulation of the kinetic 
equations of the AGN-201. Appendix B contains details of the analog 
computer set-up. 

The computer used weg a Donner Model 3100, containing 30 amplifiers 
[1] . A Donner Model 3735 electronic multiplier was used to obtain the 
product of the two time-varying quantities mentioned in Section 2, The 
3100 has 15 amplifiers suitable for use as integrators or as summer- 
inverters and 15 amplifiers which can be used as summer-inverters only. 
40 ten-turn helipot ootentionstera four of which are ungrounded, are in- 
cluded in the computer. 

‘The solution used in this investigation required 11 integrators, 15 
summer-inverters and 39 potentiometers. sua Schematic Drawing in Appen- 
dix B. 

Two programming problems, relating to magnitude scaling, made the 
A computer set-up rather intricate and unwieldly. First of all, 
the use of an inherently drift-prone electronic multiplier required 
special effort to arrange scaling of the problem so that the input volt- 
ages to the multiplier were as near as practicable to the maximum allow- 
able of 100V. This technique, including the division of the multiplier 
output by a large factor before using it in the rest of the solution, 
eliminated the problem of multiplier drift from practical consideration. 
Secondly, the wide differences between the coefficients of the different- 
ial equations required some very small potentiometer settings. In order 
to achieve maximum accuracy of simulation, additional fractional inverters 


(amplifiers 18 and 19, for example) were used to avoid potentiometer 
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settings near the extreme limits of arm travel. 
The scaling requirements NEE resulted in a limitation 
of the time duration for a computer run. When the analog was operated 
by itself without any controlling action being applied to it, overload- 
ina of Several amplifiers occurred after 20 to 30 seconds, depending upon 
the magnitude of the step disturbance used. This was not a serious limita- 
tion since the solution was not permitted to stray far from equilibrium 
when the analog was placed ín closed-loop operation with the Controllers. 
The computer equations used were derived from the kinetic equations 
previously developed. Reference 3 was used as a guide in this work. From 


Section 3, the kinetic equations are: 


of pikua- + ean] a 
+ pe” E: MOK CA 
d i ' | ` ' 
zr = Bissen — Ài C4 Ps 
i = 


= 1,2,...6 
At the steady-state equilibrium condition, i 


dů - E4913 Ci = Cilo), h= nle) 


Equation (4) then becomes, 


A y ze) -MU 


(4) 


Thus, 


nues cba 
Cilo) - [ip e.” no) (5) 


The coefficients in eq. (2) and (4) differ widely in magnitude, thus it 


is convenient to redefine the variables as follows, 
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hn C£) 


n, c h (o) 
NL MF. (6) 
Ca, T Calo) 








ok, = = ok = n (dollar units) 


Substituting (5) and (6) into (2) um (4) yields 


dn. = dei (1-8) - nys WAA CL Q^) 


= e ( $t + 1) A n, — ACA (4") 


at 
The normalized n and C can be written in the following forms, 


h = / * An, | 
( Foot note 2) (7) 
G, = / + AC; | 


By using this form, it is possible to operate only on the non-equilib- 


ium portion of the variables, This further reduces the effect of any 


Multiplier drift. Note that 
n- dan. 
dt dt 
dů: = de 


(219) 
+ 4 OB: ACK, 


dáte = M Stach &£An,$k, 上 AiAn 4") 
— M ACi, 


lye are using here the approximation that Mk H Sk, is equal 
to (k__.-1). The precise definition is (k ker The maximum error 
caused by this approximation is less than ut when Ê s = 1.0064, at prompt 
critical, 

, 2An, and AG, are variational changes of the normalized guantities 
from unity. 
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Machine voltage representation of the varíables is computed by using 


conventional scaling factors. The prime denotes machine voltages. 


L = at 
n = bh, 
CR ees 
Sk, = e Sk; 


The final form of the kinetic equations suitable for direct simulation 


on an electronic analog computer is 


Tele Tue [ratore 
- [y Jan! + [27 ES y; AC, 
ut Jan sk (9) 


Gala - Tëlee 


/ 
dali a Ai 0 g 
at 7 











The following scaling factors were found convenient for this applica- 


tion, 


1 volt/sec. 


a = 
b - 10 volts/normalized neutro: density. 
g. =g = 100 volts/normalized precursor density 


i 
= 20 volts/dollar reactivity. 


; / Ao 
Because the An, and AC, form of the variables was used, the initial 


onditions were zero for all the integrators. 


The conventional transfer function for a motor driven device ís 


Where p represents the time constant of the device. Since 


| 
Size)? 


equipment was not available for measuring E, and since the rods on the 
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AGN are small and are driven through a high gear ratio, an estimated 
value of 0.1 was assigned to E . It was known that the control rod 
was operating on the linear portion of the rod calibration (reactivity 


vs rod position) curve, See Appendix B for details of the rod simula- 


tíon, Using the measured speed of travel and knowing the linear worth 


of the rods, the reactivity rate of the two control rods was determined. 


Fine Control Rod bk /sec. 
Slow Speed 0.282 x 104 
Fast Speed 0.860 x 107^ 
Coarse Control Rod 
Slow Speed 1.71 x 10^ 
Fast Speed | 2.90 x 10^ 


The computer scaled quantity becomes 
e e 4 一 
TE Y ocn = od x Io" Sl. 7 Ce 
Two checks were made on the accuracy of simulation, Since the 
design of the control system was based on the performance of the system 
with the analog simulation, every effort was made to achieve a simulation 
that accurately duplicated the reactor behavior. One check consisted of 
comparison of the analog solution with a digital computer solution of the 
time response of a nuclear reactor, The other check consisted of fitting 
' the analog solution to the actual measured behavior of the iu it- 
self. A description of the latter check procedure is given later in this. 
section. | | 
With regard to the digital computer check method, Cohn and Toppel 
of the Argonne National Laboratory devised a code for the IBM-704 which 
calculated reactor behavior following a step change in reactivity [5] o 


With some minor modifications, the code was found suitable for the 
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Control Data Corporation 1604 Computer. 
The code is a FORTRAN ón which solves the Inhour Equation 
for its roots by an iterative process. A LaPlace Transformation approach 
was used in setting up the solution. The Inhour Equation roots were 
then used to obtain a time solution for normalized neutron population, 


by means of the relationship; 7 wt 


h(t) _ i p 
a = E (A Bj)C + c 


where A, and B, depend upon the magnitude of the step reactivity dis- 


turbance, the inhour roots,W,,and the magnitude of the extraneous neu- 
tron source. They are E solution for ease of computation. 
C depends upon the reactivity step and the neutron source only. The time 
intervals need not be regular and were left to the choice of the program- 
mer, See Appendix A for a discussion of the manner in which the delayed 
neutron effectiveness was included in this calculation, 

The code presumes equilibrium of the delayed neutron emitters be- 
fore the introduction of the E'R step. K sedona code, devised by 
the same authors, permits consideration of the case where equilibrium has 
not been achieved, 

As mentioned above, the program makes provisions for including the 
effect of an extraneous neutron source, The AGN-201 does have a Ra-Be 
Source for start-up purposes. When this provision is used, the effect of 
the source ís put in as an initial negative excess reactivity, With a 
Source of neutrons present in a reactor, constant neutron population 
does not indicate criticality (zero excess reactivity). On the contrary, 


the reactor is subcritical by a small amount; Kee is less than unity and 


the difficiency of neutrons is made up by the source. A convenient way 
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of expressing this lack of criticality is to assume an initial value of 
negative excess reactivity. 
It was decided to exclude the source from consíderation in thís and 


other kinetic equation solutions because of the difficulty of determining 


an accurate value for its effect. To check on the error introduced by 


neglecting the external source term, an experiment was conducted. Criti- 
cal rod settings were compared for two conditions; criticality at 1 watt 
with the source in its usual location and GER at l watt with the 
source removed from the reactor, The difference in the two settings was 
so small that the experiment could not be reporduced accurately enough 
to detect it. 
Input data for the code includes; 
kiyo source effect (excluded in this analysis) 
Ro step reactivity 


* í l 
A prompt neutron lifetime (The symbols are those 


AA delayed neutron decay constants used in Reference 5) 
Ba delayed neutron fractions 


time intervals 
Figure 2 is a comparison of the analog time response with that obtain- 
ed from the digital computer. A Péactivity step disturbance of 10 x 10^ | 
reactivity units and a delayed neutron effectiveness of 1,14 were used in 
both solutions. The initial responses of the two solutions compare very 
well but differences appear at longer times. See Appendix A for an ex- 
^plaination for this behavior. 
Two modifications to the analog simulation were necessary in order 


to compare its performance with the AGN-201 reactor. The first modifica- 


tion concerned the simulation of the actual reactivity disturbance of the 
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reactor, The reactor was disturbed by inserting a rod filled with poly- 
ethylene into the glory hole. The rod was pushed in by hand, and it 

was acknowledged that this resulted in a fast ramp input of finite rise 
time rather than a true step input. If the analog simulation was to be 
compared accurately with the reactor then it was necessary to disturb 
the analog with a fast ramp in a manner duplicating the reactor input. 
This was accomplished using a ramp generator. 

The second modification was made necessary by the fact that the 
neutron detection system used in measuring the output of the reactor 
contained a noise filter consisting of a shunting capacitor and an input 
resistor to the amplifier. See Appendix C. This RC filter produced an 
additional time delay in the observation of the response of the reactor. 
The circuit components selected produced a time constant of 0.05 seconds. 
This filter was duplicated by a lag transfer function simulator. Fig. 

3 shows the —Ó added to the reactor simulation in order to make 


it comparable with the AGN-201. 


REACTOR 
RAMP FILTER OUTPUT 
—Yy ml 


ENERATOR KINETIC (MULATION 


EQUATIONS 






Fig. 3 Modification of analog computer to 
permit comparison with reactor 
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Details of Se ramp generator and the filter simulation are contaín- 
ed in Appendix B. 

For comparison purposes, a recording was made of the AGN-201 reactor 
output, following the insertion of a fast ramp-reactivity disturbance, 
Without changing control rod settings the time required for the reactor 
flux to increase by a factor of 10 was recorded. From this information, 


the reactor period was computed using the relationship between flux change 


and period, E 
$ = do E T 


With Po = 10, the period was computed from, 


E Æ 
TE feg,l0 — 23 


The reactivity was computed using the Inhour Eguation solution for the 
AGN-201. A value of 15 x 10 " was calculated. At the time this measure- 
ment was made it was recognized that this Saf of reactivity was too high, 
The period measurements were made soon after the initial fast-changing 
behavior of the reactor died out and the flux had settled down on an 
approximate exponential rise. By so doing the period was measured before 
all delayed neutron effects had disappeared. The resulting value of 
period was too low. The Inhour Equation then, gave a value of reactivity 
that was too high”. Measurements were made in this manner to avoid the 
effect of the negative temperature coefficient of the reactor, It takes 
one or two reactor periods for all the delayed neutron effects to die out 
after a step disturbance and the reactor to settle on a true asymptotic 


period. In this length of time the reactor power would have rísen high 


lŷn the Inhour Eguation, reactivity is related to reactor period in 
an inverse manner. See reference (14). 
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enough so that the negative temperature coefficient would be introducing 
negative reactivity. The amount of reactivity introduced would be 
varying with time and attempts to compensate for it in subsequent calcu- 
lations would be very difficult. It was decided to m value of 
reactivity obtained in the above manner as an initial disturbance for 
the analog computer with the understanding that it could be legitimately 
reduced to obtain close agreement buts the analog run and the actual 
reactor behavior. 

On the analog computer, the value of 15 x 10^ was used as a re- 
activity disturbance. Curves for reactivity rise times! of 0.1, 0.2, 
0.5, and 1.0 seconds were obtained. These are plotted as curves one 
through four in Figure 4. Curve five represents the measured reactor. be- 
havior. It was decided that a rise time of 0.2 seconds best matched the 
reactor at low values of time. In an effort to achieve more exact match- 
ing of behavior, the reactivity input to the analog was SS, At a 
value of 14.0 x 10^ ret was obtained. This final curve 
was so close to other curves that it was difficult to plot. Instead 
isolated points are shown as crosses. 

By comparison of reactivity obtained by period measurements with that 
obtained from critical rod settings, and from other theoretical considera- 
tions, it was estimated that the values of reactivity obtained from the 
reactor may be in error by as much as 3%. This does not account for the 


total difference shown above, It is known that A > the delayed neutron 


1Rise time as used here refers to the tíme between start of insertion 
of the reactivity disturbance and when the full value has been applied. to 
the reactor. See Appendix B. 
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effectiveness, has a sizeable effect on the behavior of a reactor at long 
times (after the initial fast rise). It is felt that this parameter, 
which was obtained unverified from Aerojet-General data, may be in error 
enough to account for much of the difference between the analog and re- 
actor behavior. 

While the above comparison was not analytically precise, it did lend 
strength to the assertion that the analog computer simulation represented 
the behavior of the AGN-201 with reasonable accuracy. Admittedly, the + 
agreement was produced by a "forced fit". The adjustment of the para- 
meters (rise time and reactivity) was logical, however, and was made in 
a direction permitted by theoretical considerations. 

Having established the accuracy of the analog simulation, it was 
then used as a tool for designing and testing the automatic control 
system, This was accomplished by connectíng the controller between ampli- 
fiers 9 and 10 as shown in Fig. 10, Appendix B. The output of amplifier 
9 was applied to the control grids of tube 1 of both the UP and DOWN 
trigger circuits in place of ys Shown in Fig. 11, Appendix C. The relays 
K24 and K25 were used to operate the relay contacts shown at the input 


of amplifier 10, 
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5. Description of the Control System. 


Figure 5 is a block diagram of the control system. A more detailed 


schematic diagram can be found in Appendix C. 
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Fig. 5 Block Diagram of AFC System 

The neutron detector produces an output proportional to the neutron 
flux of the reactor. Flux, in turn, is directly proportional to reactor 
power. The output of the detector is filtered, amplified and then fed to 
a comparison network where it is compared with a signal representing the 
reference power level. That is, the power level at which the AFC system 
is to maintain the reactor. If a difference piece between the flux and 
the reference signal, an error signal is generated. The contactor ampli- 
fier samples the error for size and polarity. If the error exceeds a 


given size the amplifier causes the rod drive mechanism to move the con- 


trol rod in the appropriate direction to reduce the error signal. Since 


25 








the control system is a discontinuous Ea an error signal may exiat at 
all times if it is small enough. As long as the error stays within a 
designated "dead band" no control rod movement takes place, The AFC 
system forms a rod controlling mechanism in parallel with the manual con- 
trols of the reactor, An override is provided so that the reactor opera- 
tor may disconnect the AFC system simply by manipulating the reactor con- 
trol console switch for either of the control rods. The — compon- 


ents are described in more detail in the following paragraphs, 


Neutron Detector and Signal Amplifier 


À gamma-ray-compensated ionization chamber filled with boron tri- 
fluoride gas was located in one of the access ports of the AGN-201. It 
was supplied with a positive voltage of 250 VDC and a negative voltage of 
160 VDC for compensation from an existing power supply of the reactor, 
Since the signal from the chamber was of the order of 107? amperes, the 
power supply cables and the cable returning the chamber signal were 
shielded throughout to reduce noise pick-up. The chamber signal was fed 
into a filter system consisting of a D. LA capacitor and a 0.5 megohm re- 
sistor (1 megohm in parallel with the 1 megohm internal input resistance 
of the amplifier). 

The components of this filter were chosen from noise considerations 
and by trial on the actual equipment. Reference [27] 1s an examination 
of the noise problems encountered in reactor instrumentation systems. 

It was felt that the lowest frequence of consistent noise would be 60 
cycles. An initial estimate of filter parameters was made by assuming 
a 40 db attenuation of 60 cycle noise. A Bode plot gave a time constant 


of 0.267 sec. (breakpoint frequency of 3.75 rad./sec.) In order to reduce 
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the time delay of signal changes, the time constant was decreased experi- 
mentally by decreasing the capacitor and resistor values with the neutron 
chamber and a recorder in operation. A satisfactory arrangement was 
achieved at a time constant of 0.05 sec. (Breakpoint frequency of 20 
rad./sec.) This value of time constant had a minor effect upon chamber 
signal changes but still accomplished adequate noise filtering. This 
value might have been decreased still further, but at 0.01 seconds too 
much noise was being passed. No further effort was made toward finding 
optimum fílter parameters, 

Amplification of the chamber signal was accomplished using a low- 


drift, high-gain amplifier, Brush Model BL-550. 


Comparison Network 


Discounting the complicated initial behavior of a reactor immediate- 
ly after a step reactivity disturbance, the reactor response, after a 
sufficiently long period of time, is basically of the exponential type 


characterized by the equation t 


— 


n=N.e T (14) 
We have seen that the reactor period, T, is related in an inverse manner 


to reactivity through the Inhour Equation, Thus, 


t 
Dez oe SL 
Ne 


For a given reactivity disturbance, n/n, is the same value, at any given 
time, regardless of the power level at which the disturbance occurred, 

In reactor control, error determination is not a matter of subtracting a 
reference value from a measured quantity. Rather it is the division of 
the measured quantity by the reference and the determinetion of the devi&- 


tion of this normalized guantity from unity, Thus, 
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n Mo + LN _ An 
No Po i ho 
A comparison network was chosen that would perform this division by no 


and would permit no to be varied according to the power level desired, 


El . Referring to Fig. 6, 
ogee ies 
37 +R 

by the principle of voltage division, provided the current to the trigger 


circuit grids is kept very small. 


TO GRIDS OF 





„26 V TRIGGER CIRCUITS 


Vi (from amplifier) 


R 


Fig. 6 Comparison Network 
it is convenient for the circuitry that follows (frigger circuits) 
to have A represent a base voltage plus an error. A base voltage of 
10 volts was chosen. 


V9 = 10 + Ve 


Vh 





Se Ve = Vy —I10 ER A0 
Vin -10(1+ FB) Is 


i+ PDR 
Let V, the reference voltage, eguêl 10(1 + PD/R), then, 


V. _ Vn — Vo as 100 Vin - Vo) 
e. E ae men p—————— 


Vo Vo 
TO 








If V. represents n, and m represents no» then the output error volt- 
age is of the required form An 
| V 10 (Vn - Vo) > = | 
e^ Vo no S 


The value of P/R was chosen to be 2. The output voltage, V , then bee 





potentiometer setting, D, becomes the variable parameter that determines 

Vy, the reference voltage representing the desired power level, 
= qo Ed p) 

This analysis presumes that Y never goes below 10 volts or above 30 volts 

(D = 1 glas v7 30 volts). This was easily arranged by the scale changes 


on the amplifier. 


Contactor Amplifier 


Two Schmitt Trigger Circuits arranged as shown in Fig. 11, Appendix 
C operate relays and perform the function of measuring the error and initi- 


ating control rod movement when required [21 , The trigger circuit, 


Se 


onditions: one 


y 


o 


shown in Fig. 7., is a network with two stable operating 
when e is below the voltage of the grid! of tube 2, the second when Y 
is above grid of tube 2. As Ye increases and approaches the voltage of 

the grid of tube 2, tube 1 is cut off, RS having been chosen so thet the 

cathode voltage is high enough to permit this action, ‘ube 2 ts above 

cutoff and is conducting. The relay is in the energized condition. As pm 


ube i begins 


ES 


A approaches a value slightly below cathode voltage, 
conduct a small amount. This reduces the voltage at the plate o3 Tl ând, 
through voltage divider action, reduces the grid voltage cof T2, This 


action starts to cut off T2 which reduces voltage across R5. As a 


i ; i , 
The term grid refers to the control grid, in this report, 
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Fig. 7 Schmitt Trigger Circuit 
result the cathode to grid voltage of Tl reduces and the cycle repeats 
itself. The action is cumulative and the trigger changes to the sad 
stable condition with Tl conducting and T2 cut off, thereby de-energizing 
the relay. R5 and R4 are chosen so that the trigger point is just above 
10 volts and so that the trigger action is as sharp as possible. RI and 
R2 are chosen by plate current considerations and R3 was determined by 
voltage division requirements. Eb was &pproximateiy 130 volts, The 


resulting effect is that when Ki rîses â small amount 


= 


O volts the 


& 
o 
6 


ove 
relay changes condition (de-energizes). This relay action can be used to 
start control rod motion in such a direction that E returns to 10 volts. 
A second circuit, operating above the trigger point can be used to detect 
négative error. In this case, as Vs goes below 10 volts tube 2 starts to 
conduct and the relay energizes, The two relays, one responding £o 
positive error, the other to negative error are used to close switches 


a 


in the power supply circuit of the control rod drive motor to move the 
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fu 


rod in the required direction, Figure 12, Appendix C shows the connec- 


tions to the existing AGN-201 motor control circuits. 
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6. Performance of the Control System, 

When installed on the AGN-201, the automatic control system proved 
to be a simple, effective device which met most of the design criteria. 
The requirement that accuracy of control be within 1% proved to be un- 
realistic for a discontinuous system of this type on this reactor. The 
difficulty was not so much in the inability of the system to respond with 
the necessary sensitivity but rather in the achievement of a proper com- 
promise between accuracy and control rod movement, When an attempt was 
made to reduce the controller dead band to 17 either side of the desired 
power level, the control rod remained in virtually continuous operation, 
The resulting wear on the mechanical parts and the heating of the motor 
by continual starting was felt to be unacceptable, By widening the dead 
band to about + 3% satisfactory operation was achieved. 

Fig. 8 shows a comparison of the analog simulation performance with 
the performance of the AGN-201 under essentially the same conditions 
(same reactivity disturbance and rod correction rate). A step (fast ramp) 
reactivity disturbance of about 10 x 104 units was used, A T 
correction rate of 0.860 x 1074 units/sec., from the fine rod, fast speed, 
was used in the reactor and was simulated on the analog computer. For 
purposes of comparison, the following two standards for measuring per- 
formance are defined; 

Delay Time - The elapsed time between insertion of 


disturbance and the time when the curve 
first returns to the dead band. 


Maximum Power Rise - The normalized rise of neutron flux. 
is, the peak value divided by the init 
value. 


That 
ial 
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Curve (a), from the analog computer, shows normalized neutron ceneit; 


“a 


po 


» & 


(n/n ). The initial value is unity. A delay time of 18 seconds and a- "a 


u 


maximum power rise of 1.23 was observed. The reactor behavior, curve (53. 
shows 15 seconds and 1.25 for these quantities. The dead band in curve 


otic- 


| 


(a) extends from line 16 to line 18. In Curve (b) it is readily 


able. The oscillations in the reactor behavior and the lack of them in 


Pay 


the analog curve points out an interesting fact observed in thís study. 
It was found that it was not possible to dupiicate the behavior of the 
analog trace within the dead band. While the analog behavicr from one rum 
to another was generally the same, small differences in the mechanical 
action of the relays in subsequent runs made the oscillations uneven 
"dnd different from run to run. In some cases (curve (a) is an example) 
the relay action was such that the trace returned to the dead band and 
remained relatively level for a long period of time before control rod 
action was called for. During other runs, the analog remained in oscilla- 
tion as the reactor did in curve (b). By the same token, in some of the 
reactor runs, the flux entered the dead band in such à manner that no 
oscillations occurred for a time. This was probably due to fortuitous 
action of the relay and some coasting of the rod after the control circuit 
was opened, with the result that the reactor was left very close to crit- 
icality when the controller action stopped. Whenever this happens it is 
by chance since the on-off action of the controiler is not axpectad to 
provide exactly the correct amount of reactivity correction. 

The small differences between the two curves, as well as ali differ 
ences observed between the analog computer and the reactor during tnis 


study can be attributed to uncertainties and inaccuracies in several 
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quantities. They are summarized below. 





1. Uncertainties in the amount of re 
in the reactor, 

2. The estimation of the control rod drive mechanism Cime 
constants. 

3. The estimation, done by forced fit, of the rise time of 


reactivity into the reactor. 





4, The use of 1.14 for Y with no precise knowledge of c 
derivation of that quantity and thus its âccuracy. 

5. Uncertainty in the value of neutron lifetime in this 
reactor. 

6. The inaccuracies inherent in any simulation. The fact 
that the behavior of a reactor is represented by a šel ní 
differential eguations is a mathematic approximation, or 
model, of a physical system. 

In the light of the above known limitations, it is felt thar ex- 
cellent agreement was achieved between predicted reactor behavior and 
the performance of the reactor. The design procedure employed in tais 


investigation appears to be an accurate and practical experime 


method of designing control systems. 


Ph 
€t 
p 


Fig. 9 curves were obtained from a recorder attached to one oi 
neutron channels of the AGN-201. Fig. 9{a) shows the behavior of ine 
-- â P M^ 3 

reactor when a reactivity correction rate of 0,282 x 10 w&e employed 


(fine rod, slow speed). This run was made at 0.1 watts, kete tine 


typical behavior of a discontinuous control system before end after cf 





disturbance. The set power level of 0.1 watts is located os 
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scale at about 23. The dead band extends fran 22 to 24, This curve 
shows a delay time of 50 seconds and a maximum power rise of 1.47. 

Curve (b), taken at 1 watt with the same input shows better performance 
of the control system. In this run, fine rod, fast speed ins used in- 
itially, to obtain fast error correction, then the rod speed was manual- 
ly shifted to slow after about 70 seconds to reduce the oscillation of 
the reactor within the dead band. Curve (b) was the best performance of 
the control system observed during this investigation. Note that the 
dead band is about the same percentage of the reference level in the two 
cases. A scale change has occurred on the neutron channel between the 
two curves, In curve (a) the dead band is about 0.008 watts, which re- 
presents 87% of the reference, 0.1 watts, In curve (bj) the dead band is 
about 0.06 watts, which is 62 of the reference, 1 watt. The closeness 
of these two percentages indicated the successful division by n. discuss- 
ed in Section 5 in regard to the design of the comparison network, The 
slightly smaller value for curve (b) can be attributed to the reduced 
noise level at 1 watt operation. The signal in curve (a), for 0.1 watts, 


is comparatively noisy so the curve appears to have a larger dead band, 
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F ` Improvement of the Control System, 

The major difficulty experienced in the operation of the control 
system was in the adjustment of the trigger circuits employed to actuate 
the relays. Referring to Fig. 11, Appendix C, the interaction between 
R4 and R5 in achíeving fast trigger action and low hysterasts" is diffi- 
cult to define precisely. The circuit parameters were chosen largely by 
trial. While satisfactory performance was achieved and the circuits were 
stable once adjustments were made, the adjusting procedure was time con- 
suming and had to be made before each series of runs. 

It is felt that more investigation is necessary into the behavior of 
the Schmitt Trigger Circuit in order to improve its performance in an 
application such as that made in this work. No reference was found, for 
example, on the effect of an inductive loading, such as the control re- 
lays, in the plate circuit of tube 2 of the trigger. 

Further isvekcigdtisn may show the desirability of using variable 
Gett rod drive, possibly adjusted by tachometer feedback. There are 
other controls that could be added to make control action smoother and 
more precise, The addition of these devices must be weighed against cost 
and additional complexity. This is particularly true when the control 


system is for a low power, low cost research reactor such as the AGN-201, 


luysteresis refers to the fact that trigger action does not occur 
at the same value of V_ when the signal is increasing as it does when 
the signal is decreasifig. 
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APPENDIX A i 
DELAYED NEUTRON EFFECTIVENESS | 
The delayed neutrons produced in the fission of vi? are emitted 
with an average energy of 0.515 MEV [5] . The prompt neutrons, on 4 [ 
the other hand, have an average energy of about 2 MEV at emission. zr 
Due to their lower energy, the delayed neutrons are less likely to leak ` 
out of the reactor during the slowing down process. This "less liklihood 
of leaking out'" can be expressed as an increased effectiveness of the 
delayed neutrons in producing fissions. A greater percentage of the 
delayed neutrons will reach thermal energy and cause further fissions 
than will their prompt counterparts. This is the reason for including 
a parameter of this type in reactor kinetic investigations. | 
The dvétaña delayed neutron effectiveness y = E d Bi d. j 
is usually greater than unity, In small, water-moderated reactors, 
this parameter may exceed unity by as much as 257 [12] . In work of 
low precision, its value is usually taken as one. Actually, there is 
an effectiveness for each of the delayed neutron groups. The individual 
values are difficult to determine, however, so an average quantity, for 
the aggregate, is generally used. 
The method of including y in kinetic equations varies, superficially, 
in the literature. The fundamental approach, however, is the same. 
Handler [15 defines y mathematically as shown ín Section 3. To be 


consistent, he also redefines koff to account for the different behavior 


lkeepin [19] states that the parameter may be less than unity ín 
a fast Pu - metal reactor. 
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of the prompt and delayed neutrons. Gwin [22] has a similar express- 


ion for keff’ 


L 
kort = (1-8) eat cash + Ee Bi E EET, 


28h 277) o 
_.(I—6) k EST. 
( aH ++ n ko © 


E 
— | $ 0-9 = £ M E L 
le since kooz kozí 
Po=1-6 + 4 Q4 
then Ica e EUTP 


ky = /| T-4*B* P 
This method of accounting for Y is entirely equivalent to those found 
in the literature [ 6] [26] . They all lead to the same time behavior 
of thermal neutrons in a reactor. This investigator has chosen to in- 
corporate the effect of delayed neutrons in the manner shown in Section 
3. ‘tiie approach permits a logical, mathematical development of the para- 
meter y . In the digital computer solution, it was convenient to use 
a modifying factor for P as did Goertzel [ 26 because of the particular 
form of the mathematical expressions used and because the FORTRAN code 
could be most easily modified by changing the $ values. In accounting 
for delayed neutron effectiveness by the method of modifying e, some 
other changes in the kinetic equations are also required in order to 
make the equations agree mathematically with the Handler development. 
Unfortunately, without making major modifications in the code, these 
changes could not be included. This is the reason that the two curves 
of Fig. 2 show large differences at long time values. The digital solu- 


tion does not accurately include the effects of delayed neutrons, and 


4Q 








thus is assumed to be in error. The direction of the error is not known. 
The digital solution does provide an approximate comparison for the 
analog for purposes of checking the accuracy of simulation, particularly 


at low time values. 
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APPENDIX B 


ANALOG COMPUTER PARAMETERS 


` 











Mathematical 
Pot Expression Value © Pot Setting Amplifier R input Feedback 
X 
1 ^ .0124 .0124 1 1M lf 
2 Ax S a 2 
A 0305 0305 1M E: 
Az 
3 A „1114 „1114 3 1M Lat 
Ay 
4 AN . 3013 „3013 A 1M mf 
入 
5 As 1.136 E 5 0.1M lef 
A | 
6 de 3.013 .3013 6 0.1 M E 
7 bri .00374 .0374 11 1M 0.1 M 
8 UN .02480 .2480 11 1M 0.1 M 
9 ipi .02235 „2235 11 1M 0.1 M 
10 P 4 . 04490 „4490 12 1M Ol1M 
11 La „01311 „1311 12 1M 0.1 M 
12 bebe .00475 .0475 “12 1M 0.1 M 
26 x | 
入 /上 
13 ` 44-(75) ¿0124 .0124 1 1M Lef 
14 12 (73) .0305 .0305 2 1M lat 
15 222 (15) „1114 „1114 3 1M pet 
| C 
À 
16 4 (75 „3013 „3013 4 1M JME 
À 
17 ZE (5 1.136 „1136 5 0.1 M laf 
a. Sé Co 3.013 .3013 6 O.1M Lut 
19 YY 0400 0400 1 1M Lut 
ea . * pn E 
9 4,8 
20 9A? (i00) .0981 .0981 2 1M Lé 
21 34% (100) .3580 .3580 3 1M jet 
22 2100) „9685 „9685 4 1M yet 
2 4 À58 (i56 . | 
3 e Goo) 3.653 „3653 5 0.1 M yet 
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A b 


Mathematical 


Pot Expression Value Pot Setting Amplifier R input Feedbac 
24 Zich: o 9.685 .9685 6 in ` Lat | 
| n | 
25 $e (uo) 0040 .0040 1 1M WE: | 
bab (10 Ee i 
26 JAY (100) ^ .0098 .0098 2 1M je f 1 
ea b ' 
g ^30 
27 LAE (100) .0358 .0358 3 1M yet 
2 JÀ y . 
8 dav’ (iv)  .0968 0968 4 1M yet 
gÀsO e 
29 «Ser (100) .3653 .3653 5 1 M a f 
YAL O : 
30 = (100) .9685 .9685 6 1M ynf 
31 l(1—6) , 4968 . 4968 13 1M 
S V NOE 
| 
32 V (15) „1140 „1140 13 1M 
33 = ¿0497 .0497 13 1M 
34 lo La .0638 .0638 13 1M 
e 
35 = St, See pg 45 
36 E; i See pg 45 
37 十 Er See pg 47 
38 一 Er See pg 47 
v 
39 No 10 


/ 
* 104 n, is input to potentiometers 
L S 
** 09 „is input to potentiometers 


! /, 
k OM, SE is input to potentiometers 
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PARAMETERS FOR REMAINING AMPLIFIERS 
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AMPLIFIER 
- No. USE Input R in Feedback 
7 m (7) dan 0.1 M 0.144 
8 INT e, 1.0 M Lag 
9 we HW 0.5 M o. if 

0.5 M 
10 INT e. 1M Lef 
l4 IW 10 As, 1M 1M 
15 IN Ion 1M 1M 
16 IN sk, 0.1 M 1M 
57 mw $k An’ 1M 1M 
18 inv SE ën 1M 0.1 M 
19 mw ¡SEAN! 1M 0.1 M 
20 INV Pot 35 1M 1M 
21 SUM SE, 1M 1M 

SE: LN 

2 mm Sk 1M 1M 
23 mw Sk 1M 0.1 M 
24 ww Sk, 1M 0.1 M 
25 sum /0An, 1M 0.1 M 
| he 0.1 M 0.1 M 
2% nc Sk. 1M 1M - 





Ramp Generator 






e 


(Y  DIODE 


It is required that the generator input a specified amount of reactivity 


in time t seconds. 


ma ol 


t 
C 3,max = Co 


| Es E Er where Cp = = Sk, 


. The output of integrator 8 becomes, 





o 7 Cb 

ca 二 Asst 
e 

C aa = 57 


35 
Let the characteristics of the diode be such that with a bias of DN as 


shown, the diode cuts off the integrator output at voltage ee (close to 


v 


6 ). Then, 
uU MC e 
Cs May 7 GE = Ca 
and 
Qs —~ C» 
5 Co 
We know that 
€: Ls Es m Co 一 Co 
Ost t 
Given the values of e.» en and t, solution for Ags and e, can be obtained, 


giving the desired reactivity input in the given rise time interval. For 
t = 0.2 sec, and Sk, = 10 x 107%, the output of inverter amplifier 19 


has the following shape. 
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VOLTS 





— E— RISE TIME 


o 0t o2 0$ O4 O5 
TIME (SECONDS) 


Filter Simulation 
It is desired to simulate the transfer function, 


G(s) E xci. where R=0.5 M 


RECS +! C ur 
SEN ! i 
G(s) 7 ,05S+]/ a 





The following analog computer set up will provide this simulation 
77 ij | 


| 
— 





and Ce Re — ,05 


do 
Let a, = 1 
a el 
o 
Re = 0,5 M 
Ry = 0.5 M 


Ce - D DAE 


Control Rod Simulation 


Simulation of the rod drive transfer function, l 
S(0.1S +1) 


accomplished by use of an integrator and a lag transfer function simula- 


tion identical to the filter simulation but with different parameters, 
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Ju 


RELAY CONTACTS 


Ry = 1M 
Rf = IM 
ape 
a = 1 
o 
= " 
Ce 7 0. luf 


Delayed Neutron and Reactor Parameters. 


The delayed neutron parameters used in this analysis were those 


due to Keepin and others [2o] o 


A,  = 0.0124 H 0.000211 

Az  = 0.0305 @, = 0.00141 

As =0.112 ` @s = 0.00126 

Ày 0.301 Qy = 0.00254 

As =1.14 Rs = 0.00740 

Ae 73.01 R, = 0.00268 
Q = 0.00643 


The reactor parameters were obtained from Aerojet-General [1] and Cooke 


DE 
A = 4,1 x 107? sec, 


Y = 1.14 
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ELECTROMC 
MULT IPL/ER 





| | S 25 Sb ew, 
FIG. LÍO SCHEMATIC OF AMALOG COMPUTER SIMULA TION 


APPENDIX C 
AUTOMATIC FLUX CONTROL SYSTEM SCHEMATICS 
Figures 11 and 12 are schematics of the trigger circuits and motor 
control circuits for the automatic flux control system. The following 
notes will aid in their understanding. 


A. Fig. 11 


1. Resistor values used ín this study were: 


DOWN and UP 
R122.2K 
R2 = 2.7 K 
R3 = 100 K 
DOWN 
R4 = 94 K (adjustable) 
R5 - 680 (ad justable) 
UP | 
R4 = 8.2 K (adjustable) 
R5 = 400 | (adjustable) 
But = 140 volts 


2. K22 and K23 were sensitive relays designed for pull-in 
at about 5 ma. The trigger circuits were found to have 
best operation with plate current at about 15 ma, The 


variable resistor shunting the relays were used to ty- 





pass the excess current. 

3. Being sensitive relays, K22 and K23 had only one set 
of contacts. They were used to control relays K24 and 
K25 which had the necessary number of contacts to 


perform switching in the motor control circuit. 
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4. The diodes around K24 and K25 were found necessary to 
reduce arcing when these relay coil circuits were 
opened. 

5. The following adjustment procedure, used for a parti- 
cular set of tubes, is offered for reference. The A 

voltages will be different for other tubes. 

a. Adjust E b to 136 volts. 
b, Adjust Ys to 10,0 volts with external power 


supply. 
c. DOWN 
(1) AdjustR4 until grid, tube 2 reads 10.0 
volts. 
(2) Adjust R5 until cathode reads 11.6 volts. 
d. UP | | 
(1) Adjust R4 until grid, tube 2 reads 7.7 
volts. 
(2) Adjust R5 until caede reads 10.4 
volts. 
The above TUE are with respect to - Esp? at the negative of 
the amplifier output. 
6. P stands for a 50 K potentiometer, R is a 25 K 
E D refers to the potentiometer setting 
(O£D=1)- 
B. Fig. 12. 
1. The wiring additions are those circuits showing con- 
tacts from K20, K21, K24 and K25, The other circuits 


were original circuits on the AGN-201. 
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2. - The 
a. 
b. 
Co 
3. K20 


4. K21 








existing relays have the following function 
K16 - Energized by lower limit switch, fine 
control rod , 
K17 - upper limit switch 
K18 - Energizes when a scram occurs in order to 
return all rod carriages to the down position. 
is a small SPDT 6 volts relay 


is a 6PDT telephone type relay rated at 110 volts. 











ION CHAMBER 






ADJUSTAPLE 
RESIS TOR 


Â 
E RELAY COIL 
W 






CIRCUITS 


f 6x21 


i 
t 
i 
) 
' 
l 














RELAY CONTACTS 
(NORMALLY CLOSED) 


PELAY CONTACTS 
(LORMALLY OPER) 


L SoLiD STATE 
T DIODE 


PS., POWER SUPPLY 









(TO MOTOR. GOAJTÉ O4. 
Y CIRCUITS 


FIG, 11. AUTOMATIC FLUX CONTROL SYSTEM -TRIGGER CIRCUITS 








NNOŽRFAST SPEED) OR 40% (SLOW SPEED) 


K21-C 
C K25-D 











KiB 
K25-C 
00 
Y 
N UP 
: ! P 24-A K25-B 
KIT. 

| | 

| 

oFF ON | 

; , Kë SEN WE k K20-A 
| re, ng | K2| | 
o l pe | | 
| PLOT ! | | 
K21-D E 
Esa | | 

i 1 
yee COARSE ROD «e — —|-4 
lo^ SWITCH i 
| . k 


FIG.12 AUTOMATIC FLUX CONTROL SYSTEM = 
FINE ROB MOTOR CONTROL CR CUTS 
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12. 


13. 
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